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Abstract — We present a novel and accurate approximation for 
the distribution of the sum of equally correlated Nakagami-m 
variates. Ascertaining on this result we study the performance of 
Equal Gain Combining (EGC) receivers, operating over equally 
correlating fading channels. Numerical results and simulations 
show the accuracy of the proposed approximation and the validity 
of the mathematical analysis. 



I. Introduction 

THE knowledge of the statistics of the sums of multiple 
signals's envelopes is important in the analytical perfor- 
mance evaluation, such as that of equal gain combining (EGC) 
systems. However, the evaluation of the probability distribu- 
tion function (PDF) and the cumulative distribution function 
(CDF) of these sums can be rather cumbersome even for 
the statistically independent Nakagami-m or Rayleigh fading 
channels [2] -[7]. An infinite series technique for computing 
the PDF of a sum of independent random variables (RVs) 
was derived in [2]. Applying this technique, the error rate 
performance of EGC systems under Nakagami fading was 
presented in [3], whereas, in [4], the problem was analyzed in 
frequency domain in terms of semi-analytical expressions with 
infinite integrals. Other two studies on EGC diversity in Nak- 
agami fading that use numerical integration over Gil-Palaez 
single infinite integral and Hermite quadrature over double 
finite-infinite integral are presented in [5] and [6], respectively. 
Closed form solutions for some modulation schemes are also 
obtained for dual and triple diversity under Rayleigh fading 
[2], [5]. All above mentioned works assumed independent 
fading channels. 

However, in real-life applications, fading among diversity 
branches is correlated, which renders the analytical analysis 
under correlated Nakagami fading with a particular practical 
interest. Since the joint PDF of multiple correlated fading 
branches is not known, the published results for EGC diversity 
in correlated fading channels deal primarily with the dual 
branch case [7]-[9], where error probabilities for binary and 
QAM signals over correlated Rayleigh channels are expressed 
in form of infinite series. 
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Only several papers address EGC in correlated fading 
with multiple order diversity. In [10], EGC performance was 
determined by approximating the moment generating function 
(MGF) of its output SNR, where the moments are deter- 
mined exactly for exponentially correlated Nakagami channels 
in terms of multi-fold infinite series. A completely novel 
approach for performance analysis of diversity combiners 
in equally correlated fading channels was proposed in [11], 
where the equally correlated Rayleigh fading channels are 
transformed into a set of conditionally independent Rician 
RVs. Based on this technique, the authors in [12] derive 
the moments of the EGC output SNR in equally correlated 
Nakagami channels in terms of the Lauricella hypergeometric 
function, and then uses them to evaluate the EGC performance 
measures, such as outage probability (as infinite series) and 
error probability (using Gaussian quadrature with weights and 
abscissas computed by solving sets of nonlinear equations). 

All of the above approaches yield to results that are some- 
what complex, not expressed in closed form, and require 
computation of infinite series, all of which is attributed to the 
inherent intricacy of the exact sum statistics. This intricacy 
can be circumvented by searching for suitable highly accurate 
approximations for a sum of arbitrary number of Nakagami 
RVs. Various simple and accurate approximations to the PDF 
of sum of independent Rayleigh, Rice and Nakagami RVs 
are proposed in [13]-[15], which then are used for analytical 
EGC performance evaluation. [15] uses the moment matching 
method to arrive at the required approximation. 

In this paper, we use the moment matching method to 
obtain highly accurate closed form PDF approximation for 
the sum of arbitrary number of non-identical equally correlated 
Nakagami RVs with arbitrary mean powers. We then apply this 
approximation to efficiently estimate the performance of EGC 
systems by avoiding many complex numerical calculations 
inherent for the methods in above mentioned previous works. 
Even approximate closed form expressions allow one to gain 
insight into system performance by considering, for example, 
large SNR or small SNR behaviors. 

II. An Accurate Approximation to the Sum of 
Equally Correlated Nakagami-tti Envelopes 
Let Z be a sum of L non-identical equally correlated 
Nakagami-m RVs, Z\, Z 2 ,~-, Z L , 
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The PDF of each envelope Zk, 1 < k < L is given by [1] 
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having an arbitrary second moment E[Zk] = £lk, 1 < k < L, 
same fading parameter m z (assumed to be positive integer) 
and same envelope correlation coefficient between each pair 
of RVs 
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where E[-], cov(-, •) and var(-) denote expectation, covariance 
and variance, respectively. 

We propose the unknown PDF of be approximated by the 
PDF of an equivalent RV defined as 



R = 



\ k=l 
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where Rk, 1 < k < L, represent a different set of L 
identical equally correlated Nakagami RVs with equal average 
powers, E[Rk] = £l R , equal fading parameters m R and equal 
correlation coefficient p R . Additionally, it is assumed that 

PR = pz = P- (5) 

Both the MGF and the PDF of R 2 had been determined in 
closed form as [16, Eqs. (42a) and (36)] 
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respectively, where ■, •) is the Kummer confluent hyper- 

geometric function [17, Eq. (9.210)]. The PDF of R is deter- 
mined by simple transformation of RVs, /ij(r) = 2rf R 2(r 2 ), 
which yields to 
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One now needs to determine fl R and m R so that ([8]l be an 
accurate approximation of the PDF of Z defined by (Q~|). For 
this, we apply the moment matching method by respectively 
matching the second and fourth moments of RVs Z and R: 



E[Z 2 }=E[R 2 } , 
E[Z 4 }^E[R 4 } . 



(9) 
(10) 



The second and the fourth moments of R are determined 
straightforwardly by using the MGF © and applying the 
moment theorem, i.e., 
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The second and fourth moments of Z are determined by 
applying the multinomial theorem and the results presented 
in [10, Eq. (21)], [12, Eq. (43)] and Appendix A, yielding 
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with 2-F\ (-,■;■; ') denoting the Gauss hypergeometric function 
[17, Eq. (9.100)]. Note that (O and dBJl are vahd only if m z 
is positive integer [12]. Using [18, Vol. 4, Eq. (3.35.7(4))] and 
the Lauricella transformation to assure convergence [22, pp. 
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121], (fT3T l is expressed in closed form as follows 
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to the useful branch signal. In the EGC receiver, the random 
phases of the branch signals are compensated (co-phased), 
equally weighted and then summed together to produce the 
decision variable. 

The envelope of the composite useful signal, denoted by Z, 
is given by (Q~|i, whereas the composite noise power is given by 
a EGC ~ LNq/2, resulting in the instantaneous output SNR 



i + (iv-iyp' 'i + (jv-iyp 

where F^(- ■ ■ ) denotes the Lauricella hypergeometric func- 
tion of N variables defined by [17, Eq. (9.19)]. Note that 
coefficients W(4), W(2,2), W(3, 1) and W(2,l,l) can be 
expressed in terms of the more familiar hypergeometric func- 
tions as per (B.l), (B.2), (B.4) and (B.6), respectively. 

Introducing (0 and (9) into ( fTTT i and (11), one obtains the 
needed parameters for the PDF approximation (O of Z in 
closed form as 
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where E[Z 2 ) and E[Z A ] are respectively determined from 
( fT3] > and ( TT4"1 >. Note that the fading parameter m R is typically 
calculated to a positive real number. 

A. Special Case: Sum of Identical Equally Correlated Nak- 
agami RVs 

Let the equally correlated Nakagami RVs Z k ,1 < k < 
L have same second moments E[Z%] = Qz (equipowered 
branches), same fading parameter mz (as positive integer) and 
same correlation coefficient p between each pair of RVs. In 
this case, (fl~3l and (fl4l are simplified by using (A. 6) into 

{L-l)T 2 {m z + l/2) 



E[Z 2 ] = LVL Z 



1 



m z T 2 (m z ) 

x 2 Fi (~,-~;m z ;p 



(19) 



E[Z X \ = 



^z(l-v^) 



m z 

-AL(L- 1)1^(3,1 
- L(L - 1)(L - 2)(L - 3) W(l, 1, 1, 1 



LW(4) + 3L(L- 1)W(2,2) 

■6L(L-l)(L-2)W(2, 1,1) 
(20) 



where the necessary coefficients W(ki, &2, k^, kn) are again 
calculated by ( TTBT l. The needed parameters for the PDF ap- 
proximation ([8j of Z are obtained from (jT7j and (18). 

III. Application in the performance analysis of 

EGC RECEIVERS 

We consider a typical i-branch EGC diversity receiver ex- 
posed to slow and flat Nakagami fading. The envelopes of the 
useful branch signals Zk are non-identical equally correlated 
Nakagami random processes with PDFs given by (fj), whereas 
their respective phases are i.i.d. uniform random processes. 
Each branch is also corrupted by additive white Gaussian noise 
(AWGN) with power spectral density Nq/2, which is added 
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where RVs G k = Z k /^/LN a , 1 < k < L, form a set of L 
non-identical equally correlated Nakagami RVs with E[G 2 ] = 
7fc /L, same fading parameters m z and same correlation coeff- 
icient p among the diversity branches. Note that 7^ = ilk/No 
denotes the average SNR in the fc-th branch. 

Using the results from Section II, it is now possible to 
approximate PDF and MGF of (|2TT > by and (O, respec- 
tively, with ilfi replaced by 7 = Qji/(LNo). These closed 
form approximations are then used to determine the outage 
probability and the error probability of i-branch EGC systems 
in correlated Nakagami fading with high accuracy. 

A. Outage Probability 

The closed form approximation of the outage probability of 
the EGC receiver (i.e. the CDF of ^egc) at threshold t is 
obtained by applying [18, Vol. 5, Eq. (2.1.3(1))] over © as 
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where S S •, ■) denotes the confluent hypergeometric 

function of two variables defined by [17, Eq. (9.261(2))]. 

B. Average Error Probability 

Comparing ([TJ and (01, it is obvious that the error per- 
formance of an EGC system can be approximated by the 
performance of an equivalent maximal ratio combining (MRC) 
system for which many closed form solutions exist. For 
example, [19] derives the error probabilities of L-branch MRC 
with coherent and non-coherent detection of binary signals 
in identical correlated Nakagami fading channels. Thus, the 
average bit error probabilities of the coherent BPSK system 
and non-coherent BFSK are respectively expressed as [19, Eq. 
(32)] and [19, Eq. (26)] 
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Fig. I. PDF approximation of a sum of L identical equally correlated 
Nakagami RVs for Q\ = 1 and <5 = 



and 



P 



1 



NC-BFSK 



M R 2{s) 



(24) 

where i 7 ^-; ; ; •) is the Appell hypergeometric function 
(as the special case of Lauricella Fa function of two variables) 
defined by [17, Eq. (9.180(2))]. 

IV. Illustrative Examples and Discussion 

In this Section, the proposed approximation for the sum of 
arbitrary number of non-identical equally correlated Nakagami 
channels is validated by Monte-Carlo simulations. The simu- 
lation of correlated Nakagami random signals is realized by 
using the method proposed in [21, Section VII]. 

In order to model the non-identical branch signals (i.e., 
unequal average branch powers and unequal average branch 
SNRs), we introduce exponentially decaying profile, modeled 
as 

fi fe = Q 1 exp(-8(k- 1)), l<fc<L, (25) 

where fij is the average power of branch 1 (k = 1) and 5 is 
the decaying exponent. Note that 5 = denotes the case of 
identical branch signals (i.e., equal branch powers and equal 
average branch SNRs). 

Fig. 1 illustrate the high accuracy of the proposed PDF 
approximation of RV (Q~|i for a large variety of fading scenarios. 

Figs. 2 and 3 illustrate the high accuracy of the equiva- 
lent BPSK MRC error probability (f23]i for evaluation of the 
approximated BPSK EGC error probability. 

Appendix A 

In order to determine E[Z 2 } and E[Z A ], we apply the multi- 
nomial theorem [20, Eq. (24.1.2)]. The second moment E[Z 2 ] 
can be extracted straightforwardly. The fourth moment E[Z 4 }, 
after using [4, Eq. (43)] and performing some mathematical 
manipulations, can be transformed to 
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It is obvious from O that W(3, 1) = W(1,S) and 
W(2,l,l) = W(l,2,l) = W(l,l,2), which directly yields 
to the result given by (fl4l >. 

Appendix B 

Using identities given by [20, Eqs. (13.6.9) and (22.3.9)], 

i-Fi(0, m z , — au) = 1, i-Fi(— 1, m z , — au) = l + au/m z and 
iPi(— 2, m z , — au) = 1 + 2au/m z + (au) 2 /(m z (\ + m z )), 
one directly obtains 
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where a = y^/(l - ^/p) > 0. Using [17, Eqs. (9.212 (1)) and 
(7.622 (1))], it is possible to obtain the following identity 
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resulting into 
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Fig. 2. BER of equivalent MRC BPSK system for m z = 2 and <5 = 



Fig. 3. BER of equivalent MRC BPSK system for m z = 2 and S = 0.5 



Using [17, Eq. (9.212 (3)), pp. 1023] and after some simple 
algebra, one obtains the following 
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Thus, 
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